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Using a covariant formalism, we construct a chiral kinetic theory Lorentz invariant to order 
O(h) which includes collisions. We find a new contribution to the particle number current due 
to the side jumps required by the conservation of angular momentum during collisions. We also 
find a conserved symmetric stress-energy tensor as well as the f7-function obeying Boltzmann’s H- 
theorem. We demonstrate their use by finding a general equilibrium solution and the values of the 
anomalous transport coefficients characterizing chiral vortical effect. 
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Introduction. —The role of chiral anomalies in the col¬ 
lective dynamics has attracted considerable attention re¬ 
cently. It has been known for some time [1, 2] that a 
chiral medium in a magnetic field or in rotation can re¬ 
spond by a current along the field or the rotation axis - 
the chiral magnetic or chiral vortical effects (CME and 
CVE). More recently, the effects of anomalies in medium 
have come up in several different experimental and the¬ 
oretical contexts. Charge-dependent correlations which 
may be driven by the CME [3] have been observed in 
heavy-ion collisions. The CVE in hydrodynamics has 
been discovered using gauge/gravity duality [4, 5]. Later, 
both CME and CVE were shown to be universally re¬ 
quired by the second law of thermodynamics [6]. The 
recent discovery of “3D graphene” [7, 8] and the possible 
observation [9] of the CME-induced negative magnetore¬ 
sistance [10] have opened a new experimental frontier for 
investigating physical consequences of anomalies. 

Despite the recent progress, the role of anomaly in ki¬ 
netic theory has not been completely understood. Ki¬ 
netic theory is essential for the understanding of nonequi- 
librium dynamics and is applicable when external fields 
are weak and collisions are rare, so that each particle 
moves along its classical trajectory most of the time. 
Recent literature focuses on the kinetic theory without 
collisions. It was shown [11] that anomaly is encoded in 
the momentum-space Berry curvature, and the action for 
such a motion has been derived microscopically [12, 13]. 
Although the action and the equations of motion are not 
manifestly relativistic, a hidden Lorentz invariance, in¬ 
volving nontrivial modifications of Lorentz transforma¬ 
tions, has been found up to order 0(h) [14]. Such mod¬ 
ifications lead to side jumps necessary to ensure angular 
momentum conservation in collisions. However, the cor¬ 
responding modifications to the collision term have not 
been found so far. 

In this Letter we supply this so far missing important 
piece of the theory. First, we introduce a simple covari¬ 
ant formalism allowing us to demonstrate Lorentz invari¬ 
ance in an elegant and straightforward manner. We then 
discover that the side jumps not only make the colli¬ 


sion integral nonlocal, but also require nontrivial contri¬ 
butions to the particle number, energy-momentum and 
entropy currents. We prove the validity of Boltzmann’s 
U-theorem, guaranteeing relaxation to equilibrium. We 
determine the values of the CVE transport coefficients 
from the kinetic theory. With the goal of understanding 
the CVE, we focus on the physics of collisions without 
external electromagnetic fields, which will be considered 
elsewhere. 

Spin and relativity of particle worldline.- First of all 
we need to generalize the side jump found in Ref. [14] to 
finite Lorentz transformations. Let us consider the angu¬ 
lar momentum tensor of a relativistic spinning particle, 

r^y-xy + r, (l) 

where S^ v is the spin. In relativistic classical mechanics 
the separation between orbital motion and internal ro¬ 
tation as well as the definition of the center of mass are 
ambiguous. One can shift x M by and, simultaneously, 
S by A v pP — tSPp v without changing . To define 
unambiguously the particle position one needs to im¬ 
pose a gauge-fixing condition on S^ v . For a massless 
particle (p ■ p = 0), the only Lorentz-covariant condition 
Pfl S^ = 0 is not sufficient - leaving residual shifts A 
satisfying A ■ p = 0. To fix the gauge completely one 
chooses an arbitrary frame and uses its 4-velocity n to 
impose 


= 0 , 


( 2 ) 


i.e., one requires that S^ has only spatial components 
in the frame n. Together the two conditions p,jS llv = 
n v S^ v = 0 fix S^ in terms of n and p up to an overall 
factor 


str = a 


e^ a0 Pa n 0 

p • n 


( 3 ) 


In the frame where n M = (1,0), = Ae I;,fe p fc /|p| and 

thus A is the helicity of the particle. 

The frame dependence of the spin tensor S!ff in Eq. (3) 
implies that the particle position x also depends on the 
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frame choice, so that the total angular momentum in 
Eq. (1) does not. This means that if one changes the 
frame from n to n' the position shifts, x' = x + A nn <, so 
that 

s%-s» v =rKn>-p v Kn" ( 4 ) 

Dotting this equation with n v and choosing the point on 
the shifted world line where the shift is spatial in frame n, 
A nn' • n = 0, we find 

AM S£n v e^p a npn\ 

nn ' p ■ n (p ■ n)(p ■ n') 

This is the finite generalization of the infinitesimal side 
jump found in Ref. [14]. Finite side jumps have been also 
recently considered in Refs. [15, 16]. 

Collisionless current. —We now consider kinetic the¬ 
ory, where the system is characterized by the phase space 
particle density /. As the particle positions depend on 
the frame, so will /. Let us first ignore collisions, in which 
case / is constant along the world lines. Assuming / and 
f in two frames n and n' are related by f'(x') = /( x), 
we find to linear order in ft, with A ~ O(h), 

f'(x) — f(x) = —A • df (collisionless), (6) 
where A = A. 

The naive phase-space current p^f is thus not a 
Lorentz vector, since / is not a scalar field (its value 
at a given point depends on the frame). One part of 
the solution was found in Ref. [14] : the covariant current 
must include a magnetization contribution which, in the 
classical picture, is caused by the intrinsic rotation of the 
particles. In our covariant notations 

j 1 * = p^f + S^ v d v f (collisionless), (7) 

where S^ = S£ v . Both / and S /J,V transform non- 
trivially under the frame change n —» n' according to 
Eqs. (4) and (6) but, after cancellations, 

j'r _ ju — _A^(p • df) (collisionless). (8) 

Thus the current in Eq. (7) is frame independent in col¬ 
lisionless kinetic theory where p ■ df = 0. 

Collisions will make the current in Eq. (7) frame de¬ 
pendent. To solve this problem, we have to step back 
and try to understand what contribution to the current 
we may have missed. Equation (8) hints that it is related 
to the side jump and is proportional to the collision rate. 

Collisions and jump current. —Let us look at the col¬ 
lisions more closely and, for simplicity, consider elastic 
2 —> 2 collisions. From the classical point of view, such 
collisions involve 2 incoming and 2 outgoing world lines. 
It is convenient to think of incoming particles as being an¬ 
nihilated and outgoing particles as being created in that 
process. For particles without spin, we can assume that 


all 4 annihilation/creation events happen at the same 
spacetime point x. The continuity of the particle current 
is obvious in this case. 

However, for a spinning particle, this cannot remain 
true in all frames, because that would contradict con¬ 
servation of angular momentum [14]. We assume here 
that for each given collision kinematics there is a special 
frame — the “no-jump frame” n — in which all four par¬ 
ticle worldlines converge to one spacetime point as in the 
spinless case. The natural choice for this special frame is 
the center of mass frame: n = (pa + Pb)/\/s, where pa 
and pb are the momenta of the incoming particles. To 
ensure continuity of the current in a given (lab) frame 
n ^ n we must include a “jump current” associated with 
the spacelike motion of each participant particle between 
the common collision spacetime point x and the parti¬ 
cle’s annihilation/creation point in the lab frame, x + A, 
where from Eq. (5) 


AM = a'" = A 

— 1 — L - X nn 


e^P^PaUpHry 

{p ‘ n)(p ■ n) 


( 9 ) 


This tunneling-like motion of the particle during the col¬ 
lision would be from £ + A to x if the particle is incoming, 
or the reverse if it is outgoing. Weighing by the proba¬ 
bility of the collision with each given kinematics, we are 
led to consider the current 


jn =p r f + S^d„f + f C abcd A", (10) 

BCD 

where we introduced short-hand notations for the usual 
Lorentz invariant integration over the phase space of the 
particles B , C, and D: 

da 

etc. and for the collision kernel 


Cabcd = Wcd^ab - Wab^cd , (12) 

where W is the rate of collisions with given momenta 
Pa = P, Pb, Pc and po- The signs of the two terms 
reflect the directions of the jump depending on whether 
A is incoming or outgoing. 

Let us now check Lorentz covariance of the current j M 
in Eq. (10) by considering a different frame n' , as we did 
before in Eq. (8). Comparing 4-vectors and j M we 
find this time 


/M _ J-M = - / + A • df) - A "(p • df) 

+ f Cabcd (A £ n > ~ A£„)• (13) 

BCD 


The last term can be transformed using Eqs. (5) and (4) 

n 


AM __ ^ n v _ AM _ „M ^ nn ' 


p • Tl 


p • n 


( 14 ) 
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The meaning of Eq. (14) is straightforward: the jump 
from h to n! equals the jump from n to n plus the jump 
from n to n' (A) up to a shift along the world line [the 
last term in Eq. (14)]. 

Substituting into Eq. (13) we observe that is in¬ 
dependent of the integration variables pb, etc., and thus 
can be taken outside of the integration. The remaining 
integral coincides with the collision rate 


Using Eqs. (10), (15) and the transformation of / in 
Eq. (16) we can also rewrite Eq. (19) as 


P'df= f CabcdU] 

BCD 


x 



B'C'D' 


d_ 

df 


Cab'cd' 


AA\. 

p • u J 


( 20 ) 


C(x;p=pa)= / Cabcd , (15) 

BCD 

and, since by kinetic equation p ■ df = C[f] + 0(H), this 
term cancels the A p(j> ■ df) term in Eq. (13) to order H. 

To cancel the last term in Eq. (14) substituted into 
Eq. (13) the distribution function must transform under 
the Lorentz transformation (in addition to the shift of 
the argument by A in Eq. (6)) as 

f' - / = -A • df + / Cabcd —— ■ (16) 

J p-n 

BCD 

The additional term in Eq. (16) compared to Eq. (6) 
accounts for the colliding particles undergoing the side 
jumps. Thus we verified that the phase-space current j ,J ' 
in Eq. (10) is Lorentz covariant provided / transforms as 
Eq. (16) and Cabcd is Lorentz invariant. 

Collision kernel. —Using Eq. (10) and n ■ A = 0 we see 
that / = n ■ j jn ■ p , i.e., naturally, the time component 
of the current, j°, divided by the particle energy in the 
frame n. Since the collision probability Wab^-cd must 
be a Lorentz scalar, i.e., independent of n, the frame- 
dependent / cannot directly determine Wab^cd as in 
|AT| 2 ,//i/b(1 — fc)( 1 — /d)- Instead, we must use the 
distribution function in a frame, associated with the col¬ 
lision itself. The most natural choice is the “no-jump” 
frame h 


Now, with the n-independent distribution function in 
Eq. (17) we can write [17] 

W A b^cd[J] = ^\M(s,t )\ 2 ( 2 tt)' 1 6 4 (pa+Pb -Pc-P d) 

x/aM1-/c)(1-M (18) 

where factor 1/2! accounts for the indistinguishability of 
the outgoing particles. Using the Lorentz covariant jP 
and C, we can write a Lorentz invariant chiral kinetic 
theory with collisions 

d-j = C[f\, (19) 


where ft! is the no-jump frame of the collision AB' o 
C'D'. In the form of Eq. (20) Lorentz invariance is not 
manifest as in Eq. (19), but the collision kernel is ex¬ 
pressed solely in terms of the distribution function / in 
the lab frame. Equations. (19) and (20) are equivalent 
to linear order in h. 

Conserved currents .—Since the underlying quantum 
theory of Weyl fermions is invariant under CPT, we must 
take into account antiparticles, which also participate in 
collisions. These can be easily incorporated by consider¬ 
ing the particle charge q = ±1 as an additional discrete 
index of the distribution function f(x,p,q), indices A, 
B, etc. as composite indices A = (pa-. Qa)-, etc. and ac¬ 
companying integration over p by summation over q. CP 
invariance implies A = q |A|. The net current of q is given 
by 



and its conservation, d^Jjf = 0, follows from Eq. (19) 
and the charge conservation in a collision: fC, r; J‘ qC = 0. 

Similarly, one can show that the following covariant 
symmetric (and traceless) tensor 

^ = ( 22 ) 
q J P 

is conserved d v T^ u = 0 due to the energy-momentum 
and angular momentum conservation in the collisions. 

Entropy current and H-theorem .— An important prop¬ 
erty of kinetic theory is the existence of the entropy—a 
functional of / which does not decrease with time. This 
is known as the LZ-theoreni, which guarantees that the 
system relaxes to equilibrium. To prove the H -theorem 
we need to find the corresponding covariant current HP 
whose divergence is non-negative. 

First let us generalize current j 1 ' to a current describing 
advection of a generic, for now, quantity H which is a 
function of the distribution function /. Following the 
same steps as for the current jP in Eq. (10) we find that 
the following current 

HP =pPH + SP^d v n + J Cabcd AP™ (23) 

BCD 


where jP is given by Eq. (10) and C is given by does not depend on the choice of the frame n to linear 
Eqs. (15), (12), (18) with / from Eq. (17). order in h. 
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Using Eq. (20) one can also show that 

d^W 1 = J Cabcd[I] -qJ ■ 

BCD 


Cabcd [/eq] = 0). Also expressing S in Eq. (32) in terms 
of S using Eq. (4) we can then write for / eq : 

(24) i 

g(f eq )=p-U+-S^n^-qY, (33) 


Furthermore, using the AB <->• CD symmetry of the 
amplitude \M\ in Eq. (18) we can write Eq. (12) as 

Cabcd = Wab^cd if ~ 1). (25) 

where 

Wcd^,ab _ fcfpj 1 ~~ Ia)( 1 ~ Ib) _ 

Wab^cd ~ IaIb( 1 - fc)( 1 - Id) ~ r ' [ ’ 

and express the current H M = f T~L M as 

= J Wab^cd (r - 1) ^ ■ ( 27 ) 

ABCD 

Now, choosing % so that dH/df = ln[(l — /)//], i.e., 
n = f lny + (l-/)lnj4y, (28) 


where U a = U a + (fl a p ~ 9 a Up))A a , fl = fl and Y = 

Y — (A-d)Y’. The dependence on the collision kinematics 
via vector h (in A according to Eq. (9)) drops out, as it 
must, if d a Up = fl a p and Y = const, which also means 
U = U and Y = Y. The distribution in Eq. (33) describes 
a rotating (shear-free) fluid. It is easy to check that / eq 
given by Eq. (33) solves kinetic equation (19). In the 
conventional notations U = (3u , where /3 = \JU • U and 

Y = /?/r. 

Chiral vortical effect. —Now, for the rotating distribu¬ 
tion in Eq. (33), we can calculate the number current 
Jf in Eq. (21) using Eq. (10). It is convenient to ex¬ 
press the distribution in the local comoving frame, i.e., 
choose n = u (we have not relied on n being coordinate 
independent). To linear order in gradients we find 

= (34) 

with Ua dpu 7 , n q = J2 q f p (P ' u)qf 0 and 


and using the A B and C -n- D symmetry of \M\ in 
Eq. (18) we can write for the divergence of the entropy 
current H M in Eq. (27) 

d^H" = \ j w cd (r - 1) Inr > 0. (29) 

ABCD 

The rate of entropy production d-H vanishes when r = 1. 
Equilibruim. — Let us denote, for convenience, 

9(f) = In (30) 

i.e., f(g) = l/(exp g + 1). In terms of g = g(f), the ratio 
r in Eq. (26) is given by 

r = exp(g A +gB ~ gc ~ 9d)- (31) 

The collision kernel in Eq. (25) vanishes if r = 1 (detailed 
balance), which happens if g is a linear combination of 
quantities conserved in the collision (energy, momentum, 
angular momentum and charge), i.e., 

g(f eq )=p-U+^S af} n a p-qY, (32) 

where S = Sn is the spin tensor in the no-jump frame 
(orbital momentum is zero), q is the charge and U, Y 
and fl a p = —flp a are coefficients (possibly ^-dependent). 
The distribution function / eq in another frame n unre¬ 
lated to the collision kinematics can be obtained by trans¬ 
formation (16), according to which g = g— (A -d)g (since 


f j ( p ■ (- 


dfo 

dg 


d 

4-7T 2 


T 2 
12 ’ 


(35) 


where -df Q /dg = / 0 (1 - / 0 ), T = 1//3 and f 0 is the 
Fermi-Dirac distribution to zeroth order in gradients, i.e., 
g(fo) = P(p- u- qn). 

Similarly, for the stress energy tensor in Eq. (22) we 
find 


T^ = wu^u u - pg^ + i T (^u v + w"/), (36) 


where w and p are the usual expressions for the enthalpy 
and pressure of the Weyl gas and 



which is twice the result found in Ref. [1] due to the 
contribution of the spin coupling to vorticity. For the 
entropy current in Eq. (23) we find 


IP* = W* + (38) 


where s = /3(w — /m) and 

= 2^ t ~ = ' ( 39 ) 

One can check that these results agree with the general 
form found in Ref. [18] required by the second law of 
thermodynamics. 
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